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SUMMARY 


A curve-fitting technique is developed which is useful in the representation of 
noise sources. Noise from a given source is tabulated in a spherical coordinate 
system at a fixed radius from the origin. The noise, represented by mean-square 
pressure or by sound pressure level, is then a function of frequency, polar direc- 
tivity angle, and a number of parameters which represent the physical properties of 
the source. The dependence of the noise on the frequency and directivity variables 
is represented by bicubic splines. A least-squares bicubic spline fit is made to 
noise data from each test of the source to find values of the noise data at the 
knots, or nodes, of a frequency-angle computational grid. The array of noise values 
at these grid points is called noise coordinates. These coordinates are associated 
with the parameters which describe the state of the noise source; that is, the coor- 
dinates are functions of the test parameters. The functional dependence of the coor- 
dinates on the test parameters is represented by a Taylor's series expansion of the 
coordinates in terms of the parameters. The coefficients in the Taylor's series, 
which are partial derivatives of the coordinates with respect to the various param- 
eters, are found by a least-squares fit to data from different tests of the noise 
source. 

The method is applied to the prediction of coaxial jet noise. The jets are 
axisymmetric so that a bicubic spline may be used to represent the frequency and 
directivity dependence. The frequency variation is represented by five intervals 
with natural (zero-curvature) boundary conditions. The directivity variation is also 
represented by five intervals but with zero-slope boundary conditions. This bicubic 
spline has 36 degrees of freedom or coordinates. There are five independent param- 
eters, or dimensionless groups, which represent the flow state of a coaxial jet. 

Each of the 36 coordinates of the bicubic spline is represented by a third-degree 
Taylor's series in the 5 jet parameters. These Taylor's series each have 56 inde- 
pendent terms, including the constant for the origin of the expansion. Data from 
over 540 jet noise tests were collected to solve for the coefficients of the Taylor's 
series. This data base supplied over 540 equations in 56 unknowns for each of the 
36 coordinates. These sets of equations were solved in the least-squares sense to 
produce a seven-dimensional curve fit to the jet noise data base: two dimensions for 

the spline and five dimensions for the Taylor's series. 


INTRODUCTION 

The purpose of this paper is to develop an empirical method for source noise 
predictions. The method is general in the sense that it may be applied to any source 
for which measured noise data are available. This method is limited herein to noise 
sources with an axis of symmetry. Extending the method to sources without symmetry 
is possible with only a small conceptual change, but there will be an increased 
requirement for data, data reduction, and computation. 

This method was developed for application to aircraft noise prediction. The 
NASA Aircraft Noise Prediction Program (ref. 1) contains empirical methods for noise 
sources on conventional turbofan aircraft. The principal aircraft noise sources are 
the jets, the engine fans, and the engine combustion chambers. The secondary air- 
craft noise sources are the engine turbines and the airframe. The aircraft engine 



sources are all assumed to be axisymmetric. The airframe noise is probably not axi- 
symmetric; however, there are not enough data to define conclusively the directivity 
of this source. Since airframe noise is usually less than engine noise, it is fre- 
quently assumed that, in the far field, the noise from the entire aircraft is axi- 
symmetric. The method given herein can thus be applied either to the aircraft com- 
ponent noise sources or to the entire aircraft. 

Empirical prediction methods have three principal features: a prediction of the 

overall noise, a prediction of the directionality of the noise, and a prediction of 
the spectrum of the noise, lliese features appear in each of the aircraft component 
prediction methods given in reference 1. The general similarity of empirical pre- 
diction methods for physically different noise sources, such as an aircraft engine 
fan and a jet, suggests that a general formula can be given without referring to the 
particular physical features of the noise source. The general formula is developed 
in this paper. 

An empirical prediction method must give an accurate description of the noise 
for all possible states of the source. For a particular source, its state is 
described by variables such as speed, temperature, and pressure. These state vari- 
ables are called parameters to distinguish them from other variables, such as fre- 
quency and direction, which are associated with the noise field but not with the 
state of the source. Most noise-producing machines operate over a limited range of 
state that is not far removed from a common state called the design point. This 
limited range makes it possible to develop a general empirical method for source 
noise prediction. 

An empirical prediction method gives the noise field in terms of the source 
parameters, or source state, and tables of constants which are uniquely related to 
the noise source. For example, jet noise may be predicted with a set of noise spec- 
trum curves (functions). These curves are usually given graphically; however, for 
predictions with a digital computer, these graphs are converted to tables which are 
interpolated to give the spectrum functions (curves). Reference 1 contains tables of 
constants for the turbofan-powered aircraft noise sources. These tables often con- 
tain several thousand constants. 

The tables used in an empirical prediction method are the end product of a data- 
reduction process which begins with measured source state and noise field data. The 
data tables used in the prediction method may be regarded as alternate forms of the 
original data. Emphasizing the word "reduction" in the data-reduction process 
implies that the smallest possible tables are desired as a basis for empirical noise 
prediction. One method of achieving the goal of smaller tables is to use higher- 
order interpolation schemes. Referring again to the example of jet noise, the tables 
given in reference 1 are intended for linear interpolation. By using a higher order 
interpolation scheme, such as cubic splines, these tables may be reduced in size; 
that is, they may be converted to equivalent tables with fewer entries. The selec- 
tion of interpolation rules is the first step in either empirical prediction or data 
reduction. Empirical prediction and data reduction are inverse processes with common 
interpolation rules. 

The goal of empirical source noise prediction is to compute a noise field which 
is close to the measured data. A mathematical norm is selected to give objectivity 
and precision to the concept of closeness. This objective definition of closeness 
may be used to compare different prediction methods such as a theoretical predic- 
tion and an empirical prediction. It is used first, however, in the data-reduction 
process in computing the tables of constants. The classical norm of distance squared 
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is used herein. This norm has the advantage of leading to linear equations in the 
data-reduction process. The resulting prediction method is said to be best in the 
least-squares sense. 

The basic elements of the present empirical method may be summarized now. The 
method uses least-squares curve fits to measured data. The interpolation rules are 
cubic splines for noise field variables, such as frequency and direction, and Taylor's 
series approximations for source parameters. Taylor's series are appropriate for 
source parameters since these do not vary far from some standard source state values. 
The splines are necessary for field variables since they have a larger range of 
values . 

The method herein is illustrated with application to the turbulent mixing noise 
of static subsonic coaxial jets. The flow state of a coaxial jet is a function of 
five source parameters. These source parameters can be chosen from a large group of 
interdependent parameters. Many jet noise theories use area ratio (of the outer 
nozzle to the inner nozzle) and the velocity and density of each stream. Empirical 
formulas (refs. 1, 2, and 3) often use area ratio, velocity, and total temperature of 
each stream. A large data base is necessary to empirically define a function of five 
variables . 

The data base contains 65 single-jet noise tests from Lockheed-Georgia Company 
(work done under Air Force Contract F3361 5-73-C-2032 from 1972 to 1975), 146 tests of 
nozzles with area ratios of 0,75 and 1.2 from Pratt & Whitney Aircraft (ref, 4), and 
212 tests from the NASA Lewis Research Center (refs. 5 and 6). The Lewis tests 
included nozzle area ratios of 1,2, 1.5, 2.0, and 3.33. The Pratt & Whitney data are 
unique in their coverage of the temperatures of both streams. The Lewis data cover 
the broadest range of source parameters, including supersonic flows in both 
streams. The National Gas Turbine Establishment in the United Kingdom contributed 
three data sets designated as data sets A, B, and C. The NGTE data set C has the 
largest range of area ratio, from 1,4 to 8.1, of any data set but has an ambient 
temperature outer stream. The Societe Nationale d' Etude et de Construction de 
Motuers d' Aviation in France also contributed a subsonic jet data set. This data 
base has over 800 separate jet noise tests with around 200 000 individual sound 
pressure level (SPL) measurements. 

The National Aeronautics and Space Administration (NASA) has a continuing 
research program for the development of noise prediction methods. The first jet 
noise prediction methods developed under this program were by Stone (ref, 7) . Stone 
has produced prediction formulas for a wide range of nozzle types - including 
rectangular nozzles and thrust reversers - and flow states. Stone's formulas have 
served as the best coaxial noise predictions of NASA for over 6 years. They are a 
blend of classical theory and keen personal insight derived from a regrettably small 
amount of jet noise data. 

The inverted-flow nozzle, a concept developed for possible application to 
supersonic transport aircraft, was not covered by the 1974 prediction methods of 
Stone. Therefore, Stone developed separate formulas for this flow state (ref, 8) 
with the Lewis data base (refs. 5 and 6). In this latter method. Stone used the idea 
of two spectra, one from the outer premerged stream and one from the merged stream. 
Stone's inverted-flow method is similar to Jaeck's method (refs. 2 and 3) in the 
application of this idea of two additive spectra. About the same time, Pao (ref. 9) 
and Russell (ref, 10) developed a formula for inverted-flow jet mixing noise. Pao 
and Russell used the Pratt & Whitney data set and, like Stone and Jaeck, assumed an 
addition of two spectra. The Stone and Pao-Russell methods differ in the following 
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way: Stone defined the premerged and merged spectra through modifications of the 

empirical spectral curves from his 1974 method, whereas Pao and Russell used the 
rule that the premerged and merged spectra have similar shapes but different peak 
locations and amplitudes. All these methods are included in the NASA Aircraft Noise 
Prediction Program (ANOPP) (ref. 1) so that there is some ambiguity in the pre- 
dictions of this program for the case of inverted-flow jets. The method developed 
herein is intended to replace these present ANOPP methods for coaxial jet noise. 

This paper is divided into three major sections. The first section gives some 
preliminary definitions and notations which are helpful later. This first section 
also gives four strategies for making empirical noise predictions. Data reduction 
is discussed in the second section. It is shown how the data-reduction process is 
related to the format of the noise prediction method being developed. In the third 
section, one of these prediction strategies is applied to the data base for coaxial 
jet noise, and an explicit prediction method is derived for the data base. The third 
section concludes with a discussion of the concept of an optimum coaxial jet. 


SYMBOLS 

A area, m^ 

c speed of sound, m/s 

p 2 

Cp specific heat at constant pressure, m /K-s 

D{ 0) directivity of power 

order of derivative in Taylor's series 
CD directivity of power in frequency band 

d diameter of jet, m 

E cubic spline basis function 

f frequency, Hz 

J energy flux 

J energy flux in frequency band 


k,m,m' ,n 

L 

M 

• 

m 

N 

P 


integers 


level 

number of basis functions in spline fit on frequency 
mass -flow rate, kg/s 

number of basis functions in spline fit on direction 
power, W 
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distinct permutations of independent variables in derivatives in 
Taylor's series 

power in frequency band, W 

pressure. Pa 

mean-square acoustic pressure 
2 2 

gas constant, m /K s 
relative spectrum 
spherical radius 
spectrum of energy flux 

spectrum of energy flux in given direction 

temperature, K 

velocity, m/s 

derivative multiplier 

transformed prediction parameter 

prediction parameter 

ratio of specific heats 

noise level coordinate 

Kronecker delta function 

logarithmic frequency parameter 

cubic spline basis function 

polar directivity angle, deg 

derivative 

eigenvalue 

Helmholtz number 

densxty, kg/m 

azimuthal directivity angle, deg 


solid angle, sr 



Subscripts ; 


e equivalent 

ref reference 

t total 

1 primary stream 

2 secondary stream 

00 ambient 

Abbreviations : 

GALAC Lockheed-Georgia Company 

LeRC NASA Lewis Research Center 

NGTE National Gas Turbine Establishment (United Kingdom) 

P&W Pratt & Whitney Aircraft 

SNECMA Societe Nationale d' Etude et de Construction de Motuers d' Aviation (France) 
SPL sound pressure level 


EMPIRICAL SOURCE NOISE EQUATIONS 


Preliminary Definitions 


Noise is typically measured in terms of sound pressure level L<p^^> or in terms 
of mean-square pressure <p^>. These two variables are related by 


L<p >=10 log 


10 2 
^ref 


In the following discussion, it is necessary to switch frequently between these 
measures, so that a compact notation is desirable. Consequently, levels are usually 
given in units of bels (B) rather than in decibels (dB) , The factor of 10 in 
equation (1) is not used when the level is in bels. Local atmospheric quantities are 
used for reference constants to simplify further the data analysis. As an example, 
the sound pressure level can be given as 

L<P^> = log^o B (re pV) (2) 

P ^ 

•^OD OO 


Spherical coordinates as shown in figure 1 are used to specify the position of the 
noise sensor relative to the source, which is nominally at the origin of the 
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coordinate system. The surface area of the source A establishes a length scale 
which, together with the ambient density and the ambient speed of sound c^, 

can be used to form reference variables as required by the analysis. All data are 
assumed to be in the far field so that the pressure varies inversely with radius r. 
It is also assumed that any atmospheric absorption effects have been removed from the 
measured data. 

The acoustic intensity, like the mean-square pressure <p >, varies inversely 
with the square of radius; however, the acoustic energy flux per unit solid angle 
J is independent of radius in the far field of a nonabsorbing medium. It is 
expressed as 


J 


2 . 2 ^ 

r <p > 

P c 

“oo oo 


(3) 


The energy flux in a frequency band with Helmholtz number v = f^fA/c^ is denoted by 
^ ( v) so that the total flux is given by a summation 


J 


= ^(v) 

V 


(4) 


It is understood that the summation in equation (4) is over a contiguous set of 
bands. The acoustic power P is the integral over all directions of the energy flux 


P = J J dQ (5) 

where dQ is the differential solid angle. Equation (5) h^lds for bands as well as 
total quantities. The reference quantity for power is that the power 

level is 


LP = log 


10 



B (i 


P c; 
^00 00 




( 6 ) 


Equation (5) may be used for the energy flux level since the solid angle, measured in 
steradians, is dimensionless. 


Prediction Equations 

Empirical methods of noise prediction, which predominate in current practice, 
are essentially curve fits to measured data for various noise sources. The source 
state is characterized by a set of parameters otj^ which control the amount of noise 
generated. A parameter a could be the throttle position of an engine or the 
velocity of a jet. A prediction variable such as the band energy flux J(v) may 
depend on the direction in which the noise is emitted from the source as well as a 
large number of source parameters. 

Purely empirical prediction depends on multidimensional curve fitting. The 
curve fits may be made directly to the prediction variable; however, it is sometimes 
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useful to define related variables to separate the prediction process into several 
steps. Three methods of subdividing an empirical prediction equation for the band 
energy flux J(v)/ which is equivalent to the band mean-square pressure, are given 
herein. Ihe effects of band number and direction are considered in this section, and 
the influence of the parameters is suppressed. The effect of the source parameters 
is considered in a subsequent section. 

Prediction equations may be based on the following identities: 


Jl(v, 0) 


J(9) c^(v,9) 
P J( 0) 


(7) 


V, 9) 


P ‘P(v) ^(v>9) 

p q>( v) 


( 8 ) 


J(v, 9) 


J(9) q>(v) P J(v,Q) 

P P g>(v) J(0) 


(9) 


The factor J(9)/P 
spherical directions 9. 


is a measure of the distribution of the energy flux over the 
It is denoted by the symbol D(9); that is 


D( 0) 


J(9) 

P 


( 10 ) 


and, from its definition, satisfies the condition 


J D( 9) dQ = 1 (11) 

Q 


The directivity level LD is a negative number because of the constraint equation 
(eq. (11)). For axisymmetric sources which are assumed here, the solid angle is 


dQ = 2n sin 9 d0 


( 12 ) 


The directivity index is a comparison of the directed energy flux to the average 
energy flux through the spherical surface around the source. The factor q5(v)/P is 
the power spectrum denoted by S(v) 


S( v) 


-•P( v) 
P 


(13) 


This factor is a measure of the way the acoustic power is distributed over frequency 
bands and satisfies the condition 

^ S( v) = 1 (14) 

V 
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The factor J Q)/J(.Q) is the spectrum of the energy flux in a given direction 


S(v, 0) 


■^ ( V, 9 ) 
J(0) 


( 15 ) 


and the ratio ( v, 9) /T’( v) is the directivity of the power in a frequency band, 

^(v,9) = (16) 

qs( v) 

These variables satisfy integral conditions 

^ S(v, 9) = 1 (for all 9) (17) 

V 

and 

/ ‘D(v,9) dQ = 1 (for all v) (18) 

S3 

The most complicated group in equation (9) is a measure of the difference between the 
spectrum and directivity factors. 


93( V, 9) 


P 

9>(v) J(9) 


(19) 


By associating different variables in equation (19), it follows that 


95( V, 0) 


^(v,9) 

D(9) 


( 20 ) 


and 


95( V, 0) 


5(v,9) 
S( v) 


( 21 ) 


The form chosen in equation (20) shows that is a measure of the difference 

between the directivity of an energy flux band and the directivity of the total 
energy flux. The form in equation (21) shows that 95 is a measure of the differ- 
ence between the energy flux spectrum in a given direction and the power spectrum. 

In terms of these variables, three possible prediction equations are 


Jiv, 9) = P D(9) fi(v,9) 


Jl(v,0) = P S(v) 9)(v,9) 


( 22 ) 

(23) 
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j/(v,e) = p D(0) s(v) <i2(v,e) 


(24) 


It is a matter of individual preference which, if any, of these prediction equations 
are used in making curve fits to the data. They have the advantage of separating 
effects so that different qualities of the noise source can be studied independently. 
It is easier to consider a single integrated measure of noise, such as total power 
P, than to look collectively at J(v,6), a variable which is two-dimensional for 
fixed values of the source parameters. On the other hand, the total prediction 
formula may eventually require all dimensions, such as introduced by 
“12 (v, 6) in equation (24). 

It is possible that the number of source parameters affecting Q(v,0) may be 
less than the number affecting the power P; therefore, some savings in the number of 
dimensions are possible in the curve-fitting process. The use of the separated terms 
in the prediction equation also facilitates approximation. Many empirical formulas 
are based on the assumption that ^ = 1 . Ihe physical interpretation of this 
mathematical condition is that directivity effects and spectral effects are 
separable. 

The constraint conditions on the separate terms in the prediction equation must 
always be imposed and can be troublesome, as is seen later. The relative spectrum 
^(v,0) satisfies three constraints, not all of which are independent. 


j* D(0) <12(v,0) dQ = 1 (for all v) 


(25) 


J]s(v) Q(v,0) = 1 
V 


(for all 0) 


(26) 


D(0) R(v,0) dQ = 1 


(27) 


In terms of a prediction of the sound pressure level in decibels, the prediction 
equation (eq. (24)) becomes 


A 

SPL(v,0) = 10 log^g Jl(v,0) + 10 log^Q ~2 + 20 log^^ 

r ^ref 


dB 


[28) 


This sound pressure level is corrected easily to standard atmospheric conditions by 
using the value 197 for the last term in equation (28). The level LJ in equa- 
tion (28) is derived from equation (24) as 


Lj(v,Q) = log.,, P + log D(0) + log S(v) + log,„ 92(v,0) B (re p c A) 

lU lU 1U 'oooo 


( 29 ) 
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DATA REDUCTION 


It is assumed that the preliminary steps of removing absorption and reflection 
effects from the noise data have been completed and that the measured data are avail- 
able as energy flux spectra Mv,9) (re Po,c^A) for a number of tests of the noise 
source. The data-reduction process involves two further steps. The first step is 
to fit smooth curves through the energy flux spectra by using bicubic splines. The 
curve fits for each test are defined by a discrete set of coordinates whose values 
are established by the curve-fitting process. Each of these coordinates is poten- 
tially a function of all the source parameters. The second step in the data- 
reduction process is to fit curves, by using multidimensional Taylor's series with 
the source parameters as independent variables, to each of the coordinates defined 
by the first data-reduction step. 


Bicubic Splines 

The theory of cubic splines is well developed for one-dimensional functions. 

One method of extending this theory to multiple dimensions is to use tensor products 
of unidimensional functions. The method used here is fully developed by De Boor 
(refs. 11 and 12) for two dimensions and will be described briefly. The source is 
assumed to have an axis of symmetry so that the energy flux spectrum is a function of 
only two variables; the band Helmholtz number v and polar direction angle 0. If 
the data are given in constant-percentage bands, such as 1/3 octave, a logarithmic 
frequency variable, 

Ti = log V (30) 


is useful. 

Cubic splines are defined by De Boor in terms of basis functions. A basis func- 
tion for the frequency variable ri is defined on the interval terms of 

cubic polynomials on subintervals ( hm- i » hm ^ ' with m=1, 2, ..., M, The points 

are called nodes or Icnots. The basis functions are of class that is, they 

have continuous derivatives up to and including the second derivative. These 
continuity conditions relate the constants of the cubic polynomial for each 
subinterval so that, if the slope of the basis function is given at and the 

basis function is defined completely by these slopes and its values at the m + 1 
node points. A set of M + 1 basis functions is defined by 




( 31 ) 


together with the end slope conditions. De Boor uses the end conditions of zero 
slope for each basis function E^(ti), with m = 0, 1, . . . , M, and adds two basis 
functions: which has a unit slope at H = tIq ^M+2 '^hich has a unit 

slope at T] = Both these additional basis functions have values of zero at each 

node. A slightly different definition of the basis functions is used herein. Each 
basis function is required to satisfy the same homogeneous boundary condition at each 
end of the interval (t1o''0m^* boundary condition may involve only the derivatives 

of the function but not the function itself. The conditions of zero curvature were 
used at the ends of the range of the frequency variable, and conditions of zero slope 
were used at the ends of the range of the polar angle variable. When these boundary 
conditions are applied to each basis function, there is the same number M + 1 of 


basis functions as there are node points. Figure 2 depicts basis functions 
as defined herein. 


Noise Level Coordinates 


Since the noise variables such as <p > and J are positive, the noise vari- 
ables such as S(n) must be positive also. This condition is guaranteed if the 
curve fit is made with a noise level variable such as L<p^>, L<^, and LS. The 

following expressions for the levels of equation (29) are used: 


M 

LS(„) - E EJ„) 
m=0 


N 


LD(9) = V 0 (6) LD 
n n 


n=0 


M N 

L'12(ti, 9) = T" V E (t)) 0 (9) Lq? 

t m n n 

m=0 n=0 


mn 


M N 

L^(n, 6) = Z Z E (rt) 0 (9) LJ 
m=0 n=0 


mn 


(32) 

(33) 

(34) 

(35) 


The level coordinates such as are the levels at the node points of the splines 

as shown in figure 3. Assuming that L^(ti. ,9.), with i = 1, 2,..., I>M and 
j = 1, 2,..., J>N, has been measured in a tesi, equation (35) can be solved in the 
least-squares sense for the band energy flux levels standardized node 

points total power level is calculated by 


LP = r ^ anti log [l J ( ti,9) ] dQ (36) 

J V 

Q 

where 1>3 is given by equation (35). The power level at each band node point 
is 


J* antilog[L<) (t)^,9) ] dQ (37) 

Q 

and the energy flux at each direction node point 9^^ is 

LJ = log^ Q ^antilog[L J(t], 9^) ] (38) 

V 

where 1 ^( 11 ^^^, 9) and L<^(T],9yj) are given by simplified forms of equation (35) which 
hold at the node points. 
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The node point levels for each of the variables in the prediction equations are 
then computed by simple algebraic formulas as follows: 


LD„ = LJ - LP 
n n 


LS = LP - LP 
m m 


mn ‘^mn 


IfD = Lj - LS - LP 
mn mn m 


- LD„ - LS„ - LP 
mn mn n m 


(39) 

(40) 
(41 ) 

(42) 

(43) 


Coordinate Constraints 


The set of coordinates contains (M + 1)(N + 1) elements which completely 

describe the noise field for a given test through equation (35). The additional 
coordinates given by equations (39) through (43) offer alternate means of represent- 
ing the noise field by the prediction equations (22), (23), and (24); however, the 
elements of these alternate coordinates are not all independent. They must satisfy 
the integral constraints of equations (11), (14), (17), (18), (25), and (26). If 
methods were developed for predicting these alternate coordinates, there would be no 
guarantee that the predicted values would satisfy the constraint equations. An 
independent subset of each set of alternate coordinates must be selected and the 
constraint equations must be used to define the remaining elements of the set. 

Each of the summation expressions (eqs. (32) to (35)) for the level variables 
can be given as a product expression for the variables by taking anti logs of the 
equations as follows: 


M E (p) 

S(.) = II 

m=0 

N ( 9) 

_ / «« TT ^ n 



M 

N 

E^(n) 

0 (0) 


«(p, 9) 

= TT 

TT 

42 

mn 

n 

(46) 


m=0 

n=0 




M 

N 

E^(p) 

0 (0) 


<J(ti, 0) 

= TT 

TT 

m 

'^mn 

n 

(47) 


m=0 

n=0 
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The product forms should not be used for computation because repeated exponen- 
tiation and multiplication is much slower than the multiplications and additions in 
the summation forms for the levels. The product forms facilitate some proofs of 
identities which are needed in the process of satisfying the constraints on the main 
variable coordinates. Since the sum of the basis functions is a curve fit to the 
constant 1 (ref. 12, p. 110), 

M 

V E (g) = 1 (48) 

■r. 

m=0 

it follows that 

M 

TT 'V^o’ 

m=1 

The function S(g) satisfies a summation constraint (eq. (14)) so that the M + 1 
coordinates are not independent. The spectrum value Sq can be designated as 

the dependent variable with the M variables (S^^^/Sq) being independent. The value 
of Sg is found from 


M S 

LSg = -log^g ^ anti log ^ 

g m=1 0 

The independent coordinates L(Sjj^/Sg) may be used for tabulation and curve fitting 
of empirical data. Equation (50) gives the dependent coordinate LSg which in turn 
allows evaluation of all spectrum levels by the identity 

S 

LS = L -^ + LS- (m = 1,2, ...,M) (51) 

m Sq 0 

A similar procedure is used to find the directivity coordinate LDg 

-7t N 

LD = -log g 2% I sin 0 antilog ^ 0 ( 0) L — d0 (52) 

“'O n=1 " 0 

The spectrum constraints on the coordinates may be satisfied along each 

line of nodes where n is constant by using ^On the dependent coordinate and 

solving for Sg^ in terms of ^mn/^On* ^ similar procedure applies for the direc- 
tivity coordinates q) „ along each line of nodes where m is constant. 

This direct solution process is not possible for the dependent coordinates of 
the relative spectrum. However, if MN values of the array are stored, then 

the constraint equations (eqs. (25), (26), and (27)) give M + N + 1 nonlinear equa- 
tions for the missing values of the (M + 1 ) (N + 1) array Rmn* disadvantage of 

these nonlinear equations must be weighed against the advantages of the symmetric 
prediction equation (24) as compared with the prediction equations (22) and (23). 

Figure 4 shows four schemes for saving independent noise level coordinates for 
a single test. The set which is used depends on the desired prediction equation; 
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however, each set is only an alternate scheme for representing the energy flux 
spectrum level L;3(r), 9). 


Taylor's Series 

Each of the noise level coordinates is a function of the noise source param- 
eters. An empirical approximation to these level coordinate functions may be found 
by testing the noise source with a matrix of source parameters and fitting curves 
through the reduced data. A Taylor's series is a convenient multidimensional 
function for this curve fit because a noise source typically operates with source 
parameter values which are not far removed from some standard source condition. This 
standard condition is the origin of the Taylor's series expansion. Since the noise 
levels are measured on a logarithmic scale, the source parameters may also be given 
logarithmically. A factor of 100 for a source parameter range, say from 0.1 to 10, 
translates to the logarithmic range (-1,1) of this same parameter used as an argu- 
ment of the Taylor's series. 

The symbol A is used to denote any noise level coordinate. This level can 
be relative to a reference constant or can be relative to this same coordinate from 
another experiment or theory. The use of a level coordinate relative to a theoret- 
ically predicted value is a useful technique for evaluating the theory. If the 
theory correctly predicts the experimental data, then the relative coordinate A 
should have a value of 0 and the Taylor's series fit to these data will have null 
values for the constant and all derivatives. 

The independent source parameters are given by the following equation: 


a. 

1 


=‘i - "“"lO (a 

ref 1 


(53) 


so that the Taylor's series for A is 


A 




(54) 


The approximate equation (54) will be applied later to an empirical study of coaxial 
jet noise which has five independent source parameters. A third-degree approximation 
(N^ = 3) has 55 independent derivatives in the Taylor's series as shown in table I. 
The approximate series (eq. (54)) is expanded into a linear sum as 


A = 


N , 


A 

L 

x.=^ 




where 


a 

r X. X. X, 

A D • 1 


X. = 


(55) 


(56) 
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In equation (55), stands for the constant term in the Taylor's series; 

K 21 A 3 , ..., Ag stand for the first derivatives with respect to each of the five 

parameters; and is 56 to allow for the constant and all derivatives. The param- 
eter subscripts i, j, and k are functions of as given in table II. The order 

of the derivative is a count of the nonzero occurrence of the subscripts i, j, 

and k. The integer function is the number of permutations of the distinct 

subscripts. These functions are given in table II. In computing by equation 

(56), the convention is used that Xq = 1 and with i = 1, 2, ..., 5, is 

defined by equation (53). Given a number of experiments where a noise coordinate and 
the source parameters a are measured, equation (55) may be solved, usually in the 
least-squares sense, for the derivatives A„. This solution, repeated for each noise 
coordinate, gives a prediction method for the coordinates and hence for the noise 
field. 


COAXIAL JET MIXING NOISE 

A large data base has been collected to develop an empirical prediction method 
for coaxial jet mixing noise. This method is intended to replace the methods for jet 
mixing noise. The noise prediction method for a single jet is based primarily on the 
GALAC data. The noise prediction method for a coaxial jet is based on these GALAC 
data as well as the added data sets from LeRC (refs. 5 and 6 ), P&W (ref. 4), NGTE, 
and SNECMA. These data sets are briefly summarized in table III, and an indication 
is given of the amount of data used in this analysis. There are 842 jet noise tests 
in this data base; however, only the subsonic tests, 540 in all, were used herein. 


Organization of Data 

The data sets, as received, are organized in several different ways. The NGTE 
data set A has a test matrix with primary jet velocity and velocity ratio 

V 2 /V^ as independent test parameters. The primary jet total temperature T^ ^ was 
held fixed at about 700K for most tests and the secondary temperature was equal to 
ambient. Area ratios of 2, 4, and 6 were tested by NGTE. The NGTE data 

set C tests have the same independent variables with area ratios of 1 .5, 2.0, 4.0, 
and 8.0. The NGTE data set C includes inverted-flow cases where up 

to 1 .4 but the secondary stream was still at ambient temperature. 

The P&W and LeRC data are organized around nozzle pressure ratios and total 
temperatures. P&W used two area ratios, 0.75 and 1 .2; LeRC used area ratios of 1.2, 
1.5, 2.0, and 3.33. In these tests, the pressure ratio ranged from 1 to over 4.0 and 
the temperature ranged from 1 to over 3 times ambient temperature. 

The SNECMA data are organized around nozzle pressure ratio and total temperature 
with a single area ratio of 3.5. The GALAC single jet data also have nozzle pressure 
ratio and total temperature as test parameters. 

The way in which the test matrices are organized is to some degree arbitrary; 
however, there is a tendency for the data to be grouped according to area ratio, 
nozzle pressure ratio, and total temperature. The other geometric variable which 
appears is primary plug size, A^/A^ , where Aq is the plug area. The pressure and 
temperature variables are natural because they are easily monitored during experi- 
ments. The variables more often used for acoustic correlations, V/c and 
are of course functions of the experimental variables. Based on a review of the 
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available data sets, an organization of the data based on the actual experimental 
variables - and T^/T^ - appears to be more convenient. 

A review of the data base also suggests that the test variables are arranged in 
roughly geometric proportions such that they will appear evenly spaced on a logarith- 
mic scale. For example, the NGTE data set C used area ratios of 1.5, 2.0, 4.0, 
and 8.0. The number 1.5 is near The LeRC area ratios were 1.2, 1<5, 2.0, 

and 3.3. These ratios may be associated with 2^/^, 2^'^, 2 ^' , and 2 ^' . In both 

cases, the area ratios are fractional powers of 2, that is, Similar trends are 

apparent in the pressure ratio and total temperature. LeRC used pressure ratios of 
1.0, 1.2, 1.4, 1.6, 1.8, 2.2, and 3.0, with total temperature of 1, 2, 3, and 4 times 
ambient. The total pressure ratios tend to be grouped in either a subsonic range, 

0 < p. /p^ < 1 .863, or in a supersonic range, 1 .863 < p^/p^ < 3.470. There are 
usually far fewer variations of temperature than of pressure ratio. Usually only 1 
or 2 temperatures are examined at each pressure ratio. For example, the LeRC data 
contain mostly temperatures of 3 and 4 times ambient with relatively few tests having 
ambient or 2 times ambient temperature. These observations suggest that jet noise 
tests may be conveniently classified by the following procedure. 

Area ratio .- Area ratios should be assigned to a group with nominal 

va lue of 2^““^ , where n = U, 1, 2, ..., 7. The symbol U indicates undefined 
and is used to designate a single-stream jet; U can be thought of as -<» if the 
single-stream jet is considered to have an area ratio of zero. However, it is 
equally valid to regard the single-stream jet as representing all area ratios as 
long as the two streams are identical in pressure ratio and temperature. The 
geometric mean limits of each area ratio group are shown in table IV. The classi- 
fication scheme for area ratio can also be used for plug size if the plug size is 
given as Aq/A^ , where Aq is the plug area. No plug is designated by 0, a small 
plug by 1, and so on up to a high-radius-ratio case where Aq/A^ = 8. Because of the 
smaller number of tests available for this variable, it may be best to use only even- 
value indices for plug size. 

Total pressure ratio Pt/P<»*~ total pressure ratio may be grouped in rela- 

tion to the critical pressure ratio of the nozzle. If a typical value of y is 
assumed to be 1.35, the critical pressure ratio is 1.863. It is proposed to separate 
the subsonic range into four subranges which are geometrically related to the criti- 
cal pressure ratio. The boundaries of these subranges are given by fractional power 
n/4 of the critical pressure ratio and the geometric means are fractional powers 
(2n-1)/8, where n=1, 2, .... The classification scheme is shown in table V. 

Note that a pressure ratio of 1 .000 is excluded from the classification since a jet 
with pressure ratio of 1 does not exist. If a coaxial jet is operated with a pres- 
sure ratio of 1 in one stream, then it is regarded as a single jet with pressure 
ratio and temperature ratio equal to the values in the active stream. The area ratio 
is undefined and, if the active stream is the outer stream, the central nozzle is 
regarded as a plug. 

Total temperature ratio T^/Too«- Typical values of the total temperature are 1, 
2, and 3 times the ambient temperature. Some data, such as the LeRC data set, have 
temperatures to nearly 4 times ambient. In general, there are fewer temperature 
variations in the data base than in other variables; this suggests that fewer groups 
may be used in classifying the temperature. A geometrically related set of three 
groups is proposed as shown in table VI, The index indicates a geometric mean 
temperature T^/T^ = 2 ^~^ . No symbol is needed for graphic display, as is shown 
later. 
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Classification code.- The indices for geometry and stream conditions may be 
used to form a six-digit classification code for a noise test as follows: 


<vv^vy 

Geometry 


(P^ ./P„ 
t, 1 ° 


’’i- 

t/1 j 


^Pt,2/Po 


T /T ) 

t,2^ 0° ^ 


Primary 


Secondary 


The first two symbols are geometries, the next two are for the primary stream, and 
the last two represent the secondary, that is the outer, stream. An example of a 
single-stream jet code would be 


U1 3232 


This code begins with the undefined symbol U which indicates a single-stream jet 
with identical inner and outer stream conditions. The second symbol indicates a plug 
with area between 0.088 and 0.177 of the jet area. The last four digits must have 
the repeating pattern since there is only a single stream or two identical streams. 
The digit 3 indicates a total pressure ratio between 1.365 and 1.594, and the digit 2 
indicates a total temperature between 1.414 and 2.828 times ambient. 

It is possible to use the digit 0 to represent the undefined number. This may 
be more convenient for computer sorting of data. The example code would then appear 
as 


01 3232 

and is interpreted in exactly the same manner as long as it is recognized that 0 
means undefined. 

This code could be extended to include wind-tunnel or flight data by appending a 
digit for the free-stream Mach number. Let the Mach number be 1/10 of the last 
digit. Then the example jet code with a flight Mach number of 0.3 would be seven 
digits as follows: 

0132323 


Discussion of Data 

GALAC data.- The test points for the GALAC circular jet are classified graphi- 
cally in figure 5. Circular symbols are used to indicate a circular jet or a jet 
with identical streams. The number inside of the symbol shows the number of times 
that this nominal condition is repeated. Since the streams are identical, the test 
points fall on the diagonal of the figure. The majority of the 65 test points are 
in the subsonic region. Within each pressure ratio category, there are test points 
within each of the temperature categories. These data are adequate to define a 
function which is cubic in pressure ratio and quadratic in temperature. 

NGTE data sets.- The test points for NGTE data set A are shown in figure 6. 
These test points are predominantly subsonic. The outer stream is always cold and 
the inner stream always has a moderate temperature. There is a good distribution of 
area ratio shown by the multiple symbols in figure 6. The test points include 
circular jets which fall on the diagonal as well as dissimilar stream test points. 
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In general, the test matrix forms a banded pattern around the diagonal. The distance 
of the test point away from the diagonal is a measure of the difference between the 
pressure ratios of the streams. This difference variable is actually the ratio of 
the stream total pressures p^ 2^^^t 1 the coordinate system is on a geometric 

scale. The banded pattern of these ^ data may be adequate to define a function which 
is cubic in geometric mean pressure ratio and quadratic in the ratio of pressure 
ratios. The pattern of test points is repeated at area ratios of 4.0 and 6.0. These 
three area ratios may be used to define a quadratic function of area ratio. The NOTE 
data set B (fig. 7) contains test points off the diagonal which will better define 
the effect of different pressure ratios; however, only one area ratio is represented. 
The NGTE data set C (fig. 8) data have a good distribution of area ratios from about 
1.5 to 8.0. These test points define a function which is cubic in both pressure 
ratios and in area ratio. The temperatures are essentially fixed, however, with a 
moderate temperature primary stream and a cold secondary stream. 

P&W data .- The test points for the P&W data are displayed in figure 9. These 
data are unique because temperature effects were studied in detail. These points 
define quadratic variation in temperature for both streams. The four groups of tem- 
perature points at different pressure ratios may be used to define linear variations 
(in two dimensions) of the temperature effects with pressure ratio. Unfortunately, 
the effect of area ratio variation cannot be determined from these data. 

LeRC data .- The test points for the LeRC data in figure 10 show a distribution 
in all variables. The four area ratios could be used to define a cubic polynominal, 
although it may be better to limit this variable to second degree since the range of 
area ratio is only from about 1 to 3. At an area ratio of 1.2, three temperature 
test points are located at each of several pressure ratios; these three points define 
linear variations in the two temperatures. A pattern of three temperature points is 
also found at an area ratio of 3.3 so that a linear variation of the temperature 
effect with area ratio may be defined by these data. 

SNECMA data .- The test points for the SNECMA data are displayed in figure 11. 
This comparatively small data set could generate a function which is quadratic in 
geometric mean pressure ratio and linear in the temperature without coupling between 
the two effects. An independent variation in the quotient of pressure ratio could 
also be obtained from these data. 


Coaxial Jet Noise Prediction 

The form of the prediction equation selected for coaxial jet noise is obtained 
by combining equations (28) to (30) as 


SPL(t|, 9) = 10 log^Q P + 10 log^Q 0(9) + 10 log^Q S(n) + 10 log^p g2(r),9) 

7 . 


+ 10 log, „ — — + 20 log „ 

^10 2 ^10 p 

r ref 


dB 


(57) 
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The functions P, D, S, and are defined by equations (5), (10), (13), and (19) 

with the assumption of an axisymmetric noise field. The reference area A for the 
coaxial jet is defined as 

m 

A ^ (58) 

P c 

*^GO 00 

where m is the total mass-flow rate of the jet. The source to observer distance 

3 

is r and the frequency variable t) is the logarithm of the Strouhal number defined 
as 

fd 

e 

where d^ and are the diameter and velocity of the equivalent jet. 

The single equivalent jet has the same mass flow, energy flow, and thrust as the 
coaxial jet. The mass flow of the single equivalent jet is 


• • • 

m = m, + m„ (60) 

e 1 2 


where 


m = pAV 


(61) 


The conditions of equivalence of mass flow and thrust give 


V 

e 


mi Vi + m^V^ 


mi + m^ 


(62) 


Since the gas constant for air is not significantly changed by the addition of a 
small amount of combustion products, the equivalent temperature can be defined from 
the total energy flow as 


in.,[Yi/(Yi - 1)] + ~ ^ ^ 


where 




(63) 


( 64 ) 
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The equivalent ratio of specific heats is defined from the mixinq of the gases of 
each stream as 


"'^e - D] + 

Y - 1 ~ • • 

'e m^ + m^ 


(65) 


Because the fully expanded jet static pressure is equal to the ambient static 
pressure, the equivalent jet density is 



The equivalent area is 
m 

A = 
e p V 
^e e 

and the equivalent jet diameter is 



( 66 ) 


(67) 


( 68 ) 


The process for developing a coaxial jet noise prediction based on equation (57) 
for the data base is as follows: 

1. Compute the band energy f^ux in a given direction J(f,9) from the mean- 
square acoustic pressure data <p (f,9)> for each subsonic test by using 
equation (3). 

2. Convert the frequency values f to frequency parameter values n with 
equation (59) , 

3. Determine the energy flux coordinates J for standard values of n and 

mn 'iji 

9 to express J(ti,9) in the bicubic spline form of equation (35). The standard 
node point values of ti and 9 are shown in figure 3. These coordinates are 
found by solving equation (35) in a least-squares sense for the coordinates 
from the given values of g, 9, and for each test. 

4. Compute the power level LP from equation (36), the directivity coordinate 

levels LD^ from equation (39), the spectrvim coordinate levels from equa- 
tion (40) , and the relative spectrum coordinate levels f equation (43) . 

5. Apply the coordinate constraints for LS^ and as presented in 

equations (50) and (52) . 


21 



6. From the tabulated values of the nozzle pressure ratios 1 /Poo ^^1*3 

Pt 2 /P 00 ' total temperature ratios i/'^od ''^t 2 /^ 00 ' area ratio A 2 /A^ , 

build a Taylor's series expansion for each of the 36 noise level coordinates (P, 

5 values of D^, 5 values of and 25 values of series has the form 

of equation (55) , with the source parameters defined as 


^1 = 

0 

o'Pt, 

,/1 .365p ) 

1 CO 

^2 = 

log^ 

o'^t. 

2 / 1 .365p^) 

^3 = 

log^ 



^4 = 

log^ 

(T 

o' t. 



log^ 



The derivative terms to be included in the series are selected based on the 

classification codes for each independent variable. For example, it takes values 

of X. in three classification intervals with values of x„, . . . , in the same 

1 2 2 ^ 
interval to compute 5 A/3x^ . A subset of the 56 possible derivatives is selected 

in this manner. The Taylor's series expansion is then determined by solving equa- 
tion (55) with this subset of derivatives in a least-squares sense for each of the 
36 noise coordinates. 

Once the process of building the noise prediction method is completed, the fol- 
lowing steps are required to make a prediction with this method: 

1. Calculate the source parameters given by equation (69). 

2. Calculate the 36 noise coordinates using equation (55). 

3. Compute the values of S, D, and ^ for the desired values of r) and 9 

from the noise coordinates by using equations (32), (33), and (34). 

4. Compute the SPL from equation (57). 

Single jet validat ion .- This prediction method was first applied to the GALAC 
single jet data set as a test and validation of the method. Taylor's series 
expansions of the 36 noise coordinates for the subsonic cases were made by using the 
jet pressure ratio p^/p^ and total temperature ratio T^/T^. The following table 
is a summary of the standard deviations of the third-order Taylor's series fit to the 
data: 


Acoustic power, 10LP, dB 0.8 

Directivity level coordinates, lOLD^, dB 0.6 

Power spectrum coordinates, lOLS^, dB 0.6 

Relative spectrum coordinates, dB 1.5 

Average standard deviation, dB 1.1 


Based on these Taylor's series expansions, predictions were made of the SPL for 
selected tests in the GALAC data base. In addition, predictions were made with the 
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single jet method (ref. 1) for comparison. Ihe proposed prediction method should 
provide accuracy comparable to the prediction method of reference 1 for it to be a 
feasible method to apply to coaxial jets. 

A sample of the results of the single jet validation of the prediction method is 
shown in figure 12. Results are shown for directivity angles 0 of 120° and 150° 
and nominal total temperature ratios T^/T^ of 1.0 and 3.5 for various values of 
velocity ratio V/c . It can be seen from the figures that the accuracy of the two 
prediction methods is comparable. With this successful completion of the validation 
of the prediction method, the formulation of a coaxial jet prediction method was 
pursued. 

Coaxial jet validation .- The empirical source noise prediction method was 
applied to the entire subsonic coaxial jet noise data base and to each of the six 
data sets separately. First-, second-, and third-order Taylor's series expansions 
were made to compare the effect of order on prediction accuracy. A summary of the 
standard deviations for each of the coordinate fits for the Taylor's series 
expansions is given in table VII. Missing entries in the table correspond to cases 
where insufficient parameter variation exists in the data base to evaluate a series 
of that order. It can be seen from the table that the third-order fit provides a 
significant gain in accuracy over the first- and second-order fits. In addition, the 
error comparison indicates that there is a consistency of quality of fit among data 
sets . 


The prediction method based on the third-order Taylor's series expansion to all 
data sets is recommended for use at this time. The data base provided sufficient 
data to evaluate 43 terms of the Taylor's series for the 36 noise coordinates as 
given by equation (55). Table VIII gives the definition of the terms for all 
Taylor's series and table IX gives the derivative values for each of the 
36 coordinates. 

Test cases for the validation of the noise prediction method for the subsonic 
coaxial jet were selected from each data set. Comparisons were made between the 
third-order prediction based on all data sets and the third-order prediction based on 
the single data set from which the test case was selected to compare the error within 
data sets to the error among data sets. Test cases were selected to isolate the 
effects of various parameters on the noise prediction. 

The effect of primary jet temperature on the sound pressure level, with data 
from the NGTE data set A, is shown in figure 13. The data show a minor effect of 
primary jet temperature over a range of total temperature from 1.0 to 3.1. Neither 
the absolute level nor the spectral shape is significantly affected by primary jet 
total temperature. The agreement between data and predictions is good and indicates 
that these data are self-consistent and that the NGTE data set A is consistent with 
the entire data base. 

The effect of primary jet velocity on the sound pressure level, with data from 
the NGTE data set B, is shown in figure 14. The data show a significant effect of 
primary jet velocity over a range of velocity ratio from 0.7 to 1.4. The absolute 
noise level increases with primary jet velocity as expected. In addition, the 
spectral shape has a sharper peak as the primary jet velocity increases. The data 
and the individual set prediction agree quite well; this indicates a consistency with 
the data set. The difference between the prediction based on the NGTE data set B and 
the prediction based on all data sets indicates that the NGTE data set B tends to be 
lower in sound pressure level than the data base as a whole. 
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The effect of area ratio on the sound pressure level, with data from the NGTE 
data set C, is shown in figure 15. The data show a significant decrease in the sound 
pressure level as the area ratio increases from 2.0 to 8.0. This decrease is to be 
expected since the outer stream is slower than the inner stream. In addition, the 
frequency at which the peak sound pressure level occurs decreases with increasing 
area ratio; this is due to the changing equivalent jet diameter in the Strouhal 
number scaling. Agreement between the data and predictions is very good in this 
case . 


The ability of this empirical method to accurately predict noise for a wide 
range of test conditions is demonstrated in figure 16 with the LeRC data set. A 
near-circular jet is compared to two coaxial jets, one with a hot and fast inner 
stream and the other with a hot and fast outer stream. It is easily seen that the 
prediction method accurately predicts the sound pressure level over this wide range 
of test conditions. 

The effect of secondary jet temperature on the sound pressure level is shown in 
figure 17 by using data from the P&W data set. For a variation of total temperature 
ratio from 1.4 to 3.1, with other parameters held constant, little variation in the 
sound pressure level is seen. Note that the agreement between data and predictions 
is excellent in this figure. 

The effect of changing thrust on the sound pressure level is shown in figure 18 
by using the SNECMA data set. The prediction method provides accurate estimates in 
all three cases even though the noise level varies by more than 20 dB. 

The results in figures 12 to 18 show that the empirical source noise prediction 
method can be accurately and efficiently used to predict subsonic jet noise. Its 
applicability to other noise sources is yet to be determined. Furthermore, the 
Taylor's series form for the noise coordinate curve fit can be useful for further 
research into the noise problem, such as the optimum coaxial jet analysis discussed 
in the next section. 


Optimum Coaxial Jet 

Expressing noise coordinates, such as acoustic power, in a Taylor's series form 
facilitates design studies of the effect of the jet flow state parameters on the 
noise. In particular, a study of the coaxial benefit, defined herein as the 
difference in acoustic power produced by a coaxial jet and an equivalent circular 
jet, can be performed readily. The equivalent circular jet is one which has the same 
mass flow, momentum flow (thrust), and total enthalpy flow as the coaxial jet. 

A second-order Taylor's series expansion was used for the acoustic power P of 
the entire subsonic data base in terms of the equivalent jet velocity and total 

temperature T^ and the coaxial jet velocity ratio temperature ratio 

T 2 /T^ , and area ratio A 2 /A^ . There are two types of terms in this series 
expansion. The sum of terms including only or T^ or both is the acoustic 

power of the equivalent circular jet. The negative sum of all remaining terms is the 
coaxial jet benefit defined by the relation 


ALP = LP 

e 


LP = -E A X 

I' 


(70) 


24 


where the summation over X’ includes all coefficients as defined by equa- 
tion (56) which contain the parameters resulting 

expression is written in quadratic form as 

ALP = [x. ] {A . } + -^[x.] [A . .]{x.} (71) 

1 ,1 2 1 ,13 3 

where 


^1 = 


(72a) 

^2 = 

log,o(T,/2T^) 

(72b) 

^3 = 

lcq,o(V2/v,) 

(72c) 

^4 = 

1o,,„(T2/T,) 

(72d) 

^5 = 

log,o<a2/»,) 

(72e) 


The particular values of 


and A 


rll 


for this data base are 


and 



(73) 


'^ 13 ' - 


" 0 

0 

4.26 

-2.43 

-4.23" 

0 

0 

-7.16 

-36.32 

3.52 

4.26 

-7.16 

33.78 

-34.01 

19.32 

-2.43 

-36.32 

-34.01 

-15.44 

-10.34 

-4.23 

3.52 

19.32 

-10.34 

14.79_ 


(74) 


The matrices in equations (73) and (74) are partitioned to separate the equivalent 
jet terms x^ and X 2 from the coaxial jet ratio terms x^, and x^. 

Coefficients of only x^ or X 2 or both are identically zero because of the 
definition of coaxial benefit. 

Stationary values of ALP are found by setting the partial derivatives with 
respect to x^ equal to zero, that is 



(75) 
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This condition results in a set of simultaneous linear equations for the stationary 
point {x?} as follows: 


- -{4 •> 


(76) 


This stationary point may give an optimum (maximum) value for the coaxial benefit 
ALP. 


The characteristics of the stationary point of the coaxial benefit function are 
obtained by solving for the eigenvalues of the matrix [A . . ] • These eigenvalues and 
their eigenvectors are defined by ' ^ 


[A . .] [x. ] = [x. ] [\ ] 
,13 3 a 3 a a 


(77) 


[X ] = 
a 


lx. ] = 

1 a 


equation 

(74), the 

eigenvalues 

are 


-55.83 

0 

0 

0 

0 

0 

-1.34 

0 

0 

0 

0 

0 

5.50 

0 

0 

0 

0 

0 

20.84 

0 

0 

0 

0 

0 

64.05 

ors are 

r 0.01 

-0.92 

0.39 

-0.03 

0.07 

0.54 

-0.02 

-0.02 

0.81 

0.21 

0.34 

0.21 

0.34 

-0.40 

0.75 

0.77 

-0.07 

-0.1 5 

-0.40 

-0.47 

-0.01 

-0.33 

-0.84 

-0.12 

0.41 


(78) 
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The matrix of eigenvectors is orthonormal so that its transpose is equal to its 
inverse 


[x, ]'^ = [x. ] ^ (80) 

ja ja 

where the superscript T indicates transpose. The properties of the stationary 
point are found by introducing generalized coordinates {a }, which are amplitudes of 
the eigenvectors, by the transformation ^ 


{x . } = [x . ] {a } (81 ) 

3 3 « a 

This transformation allows an explicit solution for the stationary point fa^l in 
generalized coordinates as 



(4 .) 


26 


(82) 



In the present case, this solution is 



0.08 
2.79 
-1 .32 
I 0.21 

i-0.10 


( 83 ) 


and the transformation (eq. (81)) gives the stationary point in the original 
coordinates {^Cj} as 




(84) 


This stationary point is well outside the range of the present data base. The 
coordinate = -3.09 corresponds to an equivalent jet velocity of about 

10“^c^. The quadratic form for the coaxial benefit may be written in generalized 
coordinates as 


ALP 


5 

z 

a=1 


X fa 

a 




(85) 


The value of the benefit at the stationary point is the second term in equation (85) . 
If the stationary point is an optimum, or maximum value, then all eigenvalues X 
must be negative. The data of equation (78) indicate that only two of the five 
eigenvalues are negative so that the stationary point is a five-dimensional saddle 
point. Since the stationary point of the quadratic approximation is a saddle point 
outside the range of the data base, there appears to be no generally optimum coaxial 
jet. 


Figure 19 shows a contour plot of the coaxial benefit ALP as a function of 
velocity ratio and temperature ratio for nominal values of the equivalent velocity, 
equivalent total temperature, and area ratio. This contour shows a saddle point in 
the velocity ratio and temperature ratio. Contour plots at other values of equiv- 
alent velocity, equivalent temperature, and area ratio have similar saddle charac- 
teristics. The figure shows the general trend that if the velocity ratio increases, 
the coaxial benefit increases. Similarly, it shows that keeping the temperature 
ratio near 1 increases the coaxial benefit. 

Although a general search for an optimum point fails, the search technique is 
still useful for design. The jet state parameters are subject to design constraints 
which reduce the number of dimensions in the optimization problem. If certain param- 
eters are prescribed, then a stationary point may exist in the space of the remaining 
parameters. Pao (ref. 9) has estimated contours of equal coaxial benefit by plotting 
values derived from the P&W data base (ref. 4) in the space of equivalent velocity 
and velocity ratio. Pao estimated a maximum benefit of about 4 dB at an equivalent 
velocity Vg of about 1 .8c„ and a velocity ratio of V 2 /V^ near 2.0. The area 

ratio of the P&W data is near 1 and the jets were heated so that, for purposes of a 

simple comparison with Pao's result, X 2 and x^ may be set to 0 in the present 

analysis. Pao did not consider temperature ratio so that x^ is set to 0 also. The 
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quadratic form (eq. (71) ) then contains only and . A fit to the P&W data set 

using only these two parameters yields 


- {- 4 ? 70 > 


( 86 ) 


[A . .] = 

/ID 


0 

- 68.1 1 


68.11 

-73.77 


(87) 


and it is understood that i and j take on values of 1 and 3 only. 

The eiqenvalues and eigenvectors of the matrix in equation (87) are 


[X ] = 
a 


-1 1 4 .30 
0 


0 

40.57 


( 88 ) 


and 


[x, ] = 
la 


-0 .51 

0.86 


0.86 

0.51 


(89) 


The eigenvalues in equation (87) indicate a saddle point rather than a maximum. The 
generalized coordinates of this stationary point are 


- {-o.- 


0.03 


06 ^ 


(90) 


and the actual coordinates are 




(91) 


The value of the benefit at this stationary point is 0. 

The plot of the coaxial benefit contours from Pao (ref. 9) is reproduced in 
figure 20. Note that insufficient data exist to determine whether the contours 
actually close. Pao used the dashed lines to indicate that existence of an extremum 
was speculative. The location of the stationary point given in equation (91) 
corresponds to = 290 m/s and ~ ^ which is at the lower left edge of 

figure 20. A contour of the Taylor's series expansion for these data is shown in 
figure 21 . By comparing the two figures, it can be seen that the contours in fig- 
ure 20 can be interpreted to represent the positive lobe of the saddle shown in 
figure 21 . The results of both Pao and this paper clearly show that a significant 
coaxial benefit occurs for the higher values of both equivalent velocity and 

velocity ratio V 2 /V^ . This analysis shows however that there is no maximum coaxial 
benefit within the range of the data. A maximum can only be defined by the addition 
of design constraints which limit the jet state parameters to some region of the 
contour plane. Since the coaxial benefit function is a saddle, the maximum benefit 
must lie on the boundary of this design region. 
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CONCLUDING REMARKS 


An unbiased representation of a noise data base is important for empirical noise 
prediction and for validation of theoretical prediction methods. The method given in 
the present paper is one objective approach to the development of empirical 
predictions. The application of this method to the prediction of coaxial jet noise 
shows that it is a powerful and efficient technique which can represent a large 
number of physical effects. 

The combination of bicubic splines for acoustic field variables and Taylor's 
series for flow state parameters provides an economical scheme for reducing a large 
noise data base to a reasonably small table of constants. In the case of jet noise, 
a large data base with over 100 000 data elements was reduced to a table of slightly 
over 1000 constants. The data were reduced literally by a factor of 100. 

This method of data reduction and prediction requires a large data base. The 
data base must result from tests which are well distributed in the space of state 
parameters. This means that each state parameter must be independently varied in the 
test program which generates the data base. If the proper combination of state 
parameters is not available, then it will not be possible to find all derivatives in 
the Taylor's series. This was the case in the analysis of coaxial jet noise where 
only 43 of the 56 possible terms could be found out of a data base which contained 
540 tests. A carefully planned test program would have required only 56 tests, but 
the data used were gathered from independent laboratories in different nations so 
that many of the points in the parameter space were essentially repeated while others 
were overlooked. 

The methods used here are self -evaluating in the sense that the standard 
deviation of the curve fits to the data is found at the same time as the curve fit. 
The least-squares norm used in the curve fitting process is a classical measure of 
agreement which results in linear equations for the constants in the curve fits. In 
applying this method to coaxial jet noise, a standard deviation of 1 . 2 dB was found 
for the curve fit to the data base of 540 tests. 

The Taylor's series representation of state parameters facilitates optimization 
studies such as the search for a minimum noise jet given as an example herein. The 
second-order Taylor's series defines stationary points and their properties through 
classical eigenvalue analysis. Although the search for an optimum coaxial jet failed 
in the general case because the stationary point was a saddle, it still defined 
optimum directions in the state space for the reduction of coaxial jet noise. The 
minimum noise jet depends on design constraints on the jet state parameters. 

The analysis of this data set confirms the concept that inverted-flow jets are 
less noisy than equivalent single-stream jets at the higher equivalent velocities. 
Within the range of data used here, there can be as much as 6 dB benefit from the 
high-speed inverted-flow jet. On the other hand, there is a low velocity region 
where conventional coaxial jets are quieter than their equivalent single jets. The 
low speed benefit in the range of the data is limited to about 2 dB. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, YA 23665 
September 22, 1982 
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TABLE I.- INDEPENDENT DERIVATIVES OF FIVE-DIMENSIONAL 

TAYLOR'S SERIES 


First derivatives 




,2 


',3 


^4 


',5 


Second derivatives 


1 1 
21 

^ 22 
A' 

A 

31 


^,33 

41 



51 

^,52 

^,53 


,44 

,54 


,55 


Third derivatives 


A 

A' 

A' 

A' 

A' 

r 


1 1 1 
21 1 

^,221 


31 1 

‘,321 

^,331 

41 1 

‘,421 

‘,431 

51 1 

^,521 

^531 


‘,222 

A 


^,322 

^,332 


^,422 

‘,432 


,522 

^,532 


A 

A' 

A' 

t 


333 

433 

533 



,442 

^542 


,443 

^543 


,444 

^544 


A 


551 


A 


552 


A 


553 


A 


554 
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TABLE III.- COAXIAL JET NOISE DATA BASE 


Area 

Source 

ratio 

* 


U 

GALAC 
NGTE A 
NGTE B 
NGTE C 
P&W 
LeRC 
SNECMA 

0.75 

P&W 

1 .20 

P&W 

LeRC 

1 .50 

NGTE C 


LeRC 

2.00 

NGTE A 
NGTE C 
LeRC 

■ — ■ 


3.33 

LeRC 

3.52 

SNECMA 

4.00 

NGTE A 
NGTE B 
NGTE C 

6.00 

NGTE A 

8.00 

NGTE C 
Total 



* 


U indicates 


undefined. 











TABLE IV.- AREA RATIO CLASSIFICATION 


Index 


Aj/A, 


Lower limit 

Geometric mean 

Upper limit 

U 

0 

Undefined 

OO 

1 

0.088 

1/8 

0.177 

2 

0.1 77 

1/4 

0.354 

3 

0.354 

1/2 

0.707 

4 

0.707 

1 

1 .414 

5 

1 .414 

2 

2.828 

6 

2.828 

4 

5.657 

7 

5.657 

8 

1 1 .314 


TABLE V.- TOTAL PRESSURE RATIO P^/Pco CLASSIFICATION 

[y = 1 .35] 


Index 

Pt/Poo 

Lower limit 

Geometric mean 

1 

>1 .000 

1 .081 

2 

1 .168 

1 .263 

3 

1 .365 

1 .475 

4 

1 .594 

1.723 

5 

1 .863 

2.013 

6 

2.1 76 

2.352 

7 

2.542 

2.748 

8 

2.970 

3.210 


Upper limit 

1 .168 
1.365 
1 .594 
1.863 
2.176 
2.542 
2.970 
3.470 


TABLE VI.- TOTAL TEMPERATURE RATIO T./T CLASSIFICATION 

t CO 


Index 


T. /T 

t' CO 

Lower limit 

Geometric mean 

1 

0.707 

1 

2 

1 .414 

2 

3 

2.828 

4 


Upper limit 

1 .414 
2.828 
5.657 
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TABLE VII.- STANDARD DEVIATIONS OF TAYLOR'S SERIES EXPANSIONS 


Description 

All 

LeRC 

NOTE 

NOTE 

NOTE 

P&W 

SNECMA 

data sets 

data set 

data set A 

data set B 

data set C 

data set 

data set 

Third-order fit: 








Acoustic power 

1 .6 

1 .4 

1 .0 

1 .2 

1 .5 

0.7 


Overall directivity 

0.9 

1 .1 

0 .5 

0.3 

0.6 

0 .8 


Power spectrum 

1 .1 


0.5 


1 .1 

0.8 


Relative spectrum 

1 .2 


0 .5 


0.9 

1 .1 


Total 

1 .2 

1 .5 

0.5 

0.5 

0.9 

1 .0 


Second-order fit: 








Acoustic power 

2,2 

1 .7 

1 .5 

1 .4 

2.4 

0.8 

0.7 

Overall directivity 

1 .0 

1 .1 

0 .6 

0 .3 

0 .7 

1 .1 

0 .2 

Power spectrum 

1 .3 

1 .3 

0.6 


1 .1 

1 -0 

0.4 

Relative spectrum 

1 .4 

1 .6 

0 .5 


1 .0 

1 .3 

0 .3 

Total 

1 .3 

1 .5 

0.6 

0.6 

1 .1 

1 .2 

0.3 

1 

First-order fit: 








Acoustic power 

3.5 

2-4 

2.8 

3.1 

3.9 

1 .6 

0.8 

Overall directivity 

1 .4 

1 .3 

0 .9 

0 .5 

0 .8 

2 .0 

0 .2 

Power spectrum 

2.0 

1 .5 

0.9 


1 .6 

1 .5 

0.4 

Relative spectrum 

1 .7 

1 .7 

0 .7 


1 ,3 

1 .6 

0 .3 

Total 

1 .8 

1 .6 

0.9 

1 .3 

1 .5 

1 .7 

0.4 
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TABLE VIII.- DERIVATIVE MULTIPLIERS 


Definitions: 

Xi = 

='2 = 

"4 = 1°5,o(T^,2/(2TJ) 


Derivative 
index, t 

Derivative 

multiplier. 

1 

2 

1 

^1 

3 

^2 

4 


5 


6 

=^5 

7 

/2 

8 

X1X2 

9 

x,X3 

10 

X1X4 

1 1 

X,Xj 

12 

i 

X2X2/2 

1 

13 

X2X3 

14 1 

^ 2^4 

15 

X2X5 

16 

1 

X3X2/2 

17 

’' 3’'4 

1 8 

X3X5 

19 

X4X4/2 

20 

=' 4’'5 

21 

XgXs/2 

22 

X, x,x, /6 

23 

x^x^X2/2 

24 

x^x^X3/2 

25 

x,x,X4/2 

26 

x^x^X3/2 

27 

x,X2 X2/2 

28 

X1X2X3 


Derivative 
index, Si 

Derivative 

multiplier. 

29 

X^XjX^ 

30 

X 1 X 2 X 5 

31 

XiX3X3/2 

32 

x,XsX4 

33 

X1X3X5 

34 

X, x^x ^/2 

35 

X1X4X5 

36 

x,X5Xs/2 

37 


38 

X2X2X2/2 

39 

X2X2X4/2 

40 

^2^2^5^^ 

41 

X2X3X3/2 

42 

X2XjX^ 

43 

X2X3X5 

44 

X2X4X4/2 

45 

X2X^Xj 

46 

X2X5X5/2 

47 

X3X3X3/6 

48 

X3X3X4/2 

49 

X3X3X3/2 

50 

X3X4X4/2 

51 

X3X4X5 

52 

X3X5X3/2 

53 

’' 4 ’' 4 * 4 /® 

54 

X4X4X5/2 

55 

X4X5X5/2 

56 

X5X3X3/6 


IndeXj 


Derivatives values , for 10 iog^Q of - 



P 

D ( 30 ) 

D ( 60 ) 

D ( 90 ) 

D (1 20 ) 

D ( 150 ) 

S (-1 . 0 ) 

S (- 0 . 5 ) 

3 ( 0 . 0 ) 

1 

- 44,7 

- 19,6 

- 17.7 

- 15.1 

- 10.1 


- 14.1 

- 9.5 

- 15.8 

2 

64.9 

- 85.2 

- 72.3 

- 44.9 

- 27.7 


9.7 

4.5 

8.0 

3 

52.8 

-1 1 .6 

- 10.4 

- 9.3 

- 6.4 

HI 

8,3 

- 8.8 

- 2.4 

4 

27.9 

- 18.8 

- 19.6 

- 20.9 

- 16.5 

K9 

1.7 

5.6 

- 7.1 

5 

14.4 

- 9.0 

.- 10.2 

- 13.1 

- 13.1 

- 9.8 

2.9 

- 2.9 

-1 .8 

6 

-1 .4 

- 12.3 

- 12.0 

-1 2.6 

- 1 1.7 


- 7.2 

0.2 

5-5 

7 

-21 1 .3 

83.4 

132.3 

264.7 

214.7 


- 12.7 

-291 .7 

333.1 

8 

- 409.9 

-1 34.8 

- 114.7 

-81 .9 

- 13.5 


106.2 

155.5 

8,3 


135.8 

- 1 36.9 

- 96.2 

- 13,2 

35.9 


- 40.6 

- 48.5 

1 10.0 


- 80,9 

- 94.2 

- 62.6 

1 3.0 

1 1 .9 

1 3.4 

63.5 

- 3.2 

45.0 


-43 .1 

44,4 

25.5 

3,2 

7.5 

-23 .2 

-1 2 .2 

- 13.2 

13.1 


134.6 

135.5 


1 05.2 

25.1 

- 108.5 

7.6 

- 222.2 

51.8 


- 135.6 

- 2,5 


- 16,8 

- 23.8 

-1 1 .2 

32 .2 

17.9 

- 15.7 


75.3 

- 21.6 


- 3.9 

1 3.1 

- 28.8 

- 14.1 

4.4 

- 8.1 

15 

51 .4 

- 29.4 


- 7,0 

- 3.7 

- 23.0 

43 .5 

20.0 

- 15.2 

16 

23,6 

1 14.1 


15,0 

19.1 

- 35.4 

- 55.4 

-31 .0 

1 2,8 

17 

-19 ,1 

1 3,7 


- 12.6 

- 18.2 

- 17.5 

17.0 

9.8 

- 13.4 

18 

1 ,0 

- 28.5 

- 22.9 

- 10.0 

- 0.7 

- 31.3 

-1 .8 

- 1,1 

-1 .9 

19 

44,5 

- 63,9 

-63 .1 

-51 .6 

-51 .0 

10.2 

- 69.0 

- 28.3 

33 ,2 

20 

0 

0 

0 

0 

0 

0 

0 

0 

0 

21 

-1 3 .3 

- 0.8 

- 2.9 

-4 .2 

-9 .0 

- 19.9 

14,4 

-3 .0 

-7 .0 

22 

-1 326.3 

2524.2 

2825,1 

2446.1 

3722.3 

1574.2 

- 5289.4 

2337.2 

-1 1 31.7 

23 

5 145.7 

- 263.4 

-1 31 .5 

198.1 

77 .3 

- 302,6 

- 2658.9 

- 1077.2 

288.4 

24 

-2 052.9 

1370.9 

839.6 

- 354.7 

- 807.3 

- 355.7 

290.5 

1 86.0 

- 1 144,1 

25 

0 

0 

0 

0 

0 

0 

0 

0 

0 

26 

406.0 

- 23.9 

- 49.3 

- 97.4 

- 114.0 

51.6 

- 329,0 

- 87.6 

- 88.0 

27 

806.5 

518 .1 

564.1 1 

535 .4 

185.1 

40.5 

1294,9 

1190.1 

-705 .5 

28 

844.6 

- 63.7 

- 28.9 i 

40.6 

50.9 

- 65.3 1 

1 08.4 

- 96.2 

117.7 

29 

0 

0 

0 ' 

0 

0 

0 

0 

0 

0 

30 

- 388.0 

247,5 

187,9 ' 

82.2 

- 43.3 

- 30.0 

-1 5.6 

- 108.7 

- 12.7 

31 

-223 .4 

-101 .6 

- 106.8 

- 168.5 

-277 .0 

-99 .4 

233 .1 

274.6 1 

- 212.1 

32 

0 

0 

0 

0 

0 

0 

0 

0 

0 

33 

85.5 

10.0 

2.5 

- 13.9 

- 27.8 

-29 .1 

- 56.2 

-23 .1 

-65 .7 

34 

0 

0 

0 

0 

0 

0 

0 

0 

0 

35 

0 

0 

0 

0 

0 

0 

0 

0 

0 

36 

- 28,8 

- 73.5 

- 46.4 

- 17.5 

- 41.9 

1 0.5 

0 

0 

0 

37 

-15 062.4 

- 974.4 

-989 .6 

- 1102.3 

- 1502.9 

348 .0 

- 5794.5 

-1 1 53 .7 

1770,6 

38 

116.1 

1 38.1 

1 52.2 

1 61 .4 

- 18.9 

1 66,4 

1 77.7 

402.3 

- 82.5 

39 

0 

0 

0 

0 

0 

0 

0 

0 

0 

40 

609.1 

- 225.9 

- 21 2.0 

- 193,3 

- 73.0 

1 1 2.5 

- 173.8 

275.3 

- 135.0 

41 

221 .6 

- 235,0 

- 206.2 

- 184,3 

7 .6 

25 ,3 

1 7 .1 

81 .8 

-27 .4 

42 

- 193.2 

1 52.9 

122.3 

76,9 

- 33.3 

-21 .9 

111,2 

- 96.7 

64.6 

43 

-214 .3 

100,5 

69.6 

22.7 

-31 ,5 

8.9 

- 0.5 

8.4 

23.8 

44 

- 51.9 

520.5 

423.5 

274.5 

148.5 

- 244.1 

368.3 

1 68.4 

- 179.8 

45 

0 

0 

0 

0 

0 

0 

0 

0 

0 

46 

18.9 

- 1 1.5 

- 18.0 

- 27,4 

- 34,8 

17.8 

0 

0 

0 

47 

0 

0 

0 

0 

0 

0 

0 

0 

0 

48 

58.7 

64.6 

44.8 

- 0.6 

-1 .6 

- 26.9 

- 197.9 

-1 .0 

-31 .4 

49 

33.2 

- 159.1 

-111 .8 

- 34.9 

-23 .8 

2.6 

51 .4 

- 14.2 

-31 .7 

50 

18.0 

100.6 

68.6 

3.6 

32,5 

- 25.7 

-1 32.7 

7.0 

92.4 

51 

84.5 

- 79.1 

- 60.7 1 

- 33.7 

- 15.0 

-1 .3 

- 23.7 

22.6 

- 17.4 

52 

- 6.6 

- 15.4 

- 15.6 

- 23.9 

- 37,4 

40.1 

0 

0 

0 

53 

0 

0 

0 

0 

0 

0 

0 

0 

0 

54 

0 

0 

0 i 

0 

0 

0 

0 

0 

0 

55 

0 

0 

° i 

0 

0 

0 

0 

0 

0 

56 

0 

0 

0 ! 

0 

0 

0 

0 

0 

0 






























TABLE IX.- Continued 


Derivatives values, A«, for 10 log ^0 of - 


S ( 0 . 5 ) S ( 1 . 0 ) I !(- 1 . 0 , 30 ) qj (- 0 . 5 , 30 ) < U ( 0 . 0 , 30 ) ^^( 0 . 5 , 30 ) qj ( 1 . 0 , 30 ) g 5 (- 1 . 0 , 60 ) q ?(- 0 . 5 , 60 ) 


- 18.7 
-1 2.7 
26.6 
- 14.8 
0,7 
5.6 
51 7.2 
- 176.4 
84.2 

30.1 
36.9 

220.9 
- 27.5 

41.7 
- 5.2 

46.1 
- 26.4 

12.1 
0.6 
0 

- 3.7 

- 154,8 
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Fiqure 6.- Jet noise test points for NOTE data set A 
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Figure 7.- Jet noise test points for NGTE data set B 
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Figure 10.- Jet noise test points for LeRC data 
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Figure 11.- Jet noise test points for SNECMA data. 
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Figure 12,- Continued 
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Figure 14.- Effect of primary jet velocity on sound pressure level for a coaxial jet 
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e := 1 20° . 



Sound Prsssuro Level, SPL, dB Sound Pressure Level, SPL. dB 



Frequency, f, Hz 

(a) A 2 /A*| — 2«0* 


o Data 

Empirical source noise prediction 

based on all data sets 

— Empirical source noise prediction 
based on NGTE Set C only 


120 p 

u 


110 i- 



Frequency, f, Hz 

(b) A2 /At =4.3. 

Figure 15.- Effect of area 
NGTE data set C. = 

e = 120°. 


120 r- 



Frequency, f. Hz 


(c) A 2 /A 1 =8.0. 

level for a coaxial jet with 

2.8T ; T „ = 1.0T ; 

CD t, 2 ® 


ratio on the sound pressure 

1.5c ; = 0.9c ; T^ = 

00 ' 2 “ t,1 


o\ 



Sound Prissur. Leval. SPU dB Sound Preesuro L«vd. SPL <SB 


a\ 

00 


110 


100 


00 


80 



70 


I ■ I I « I I I I I I I I — I — I — I — I — I — I — I — I 

160 315 630 1250 2500 5000 10000 20000 

Frequency, f. Hz 

(a) = 0.6c^; = 0.6c^; = I.OT^; 



T ^ = 1 . OT . 
t,2 <» 


O Data 

Empirical source noise prediction 

based on ail data sets 

■ — Empirical source noise prediction 

based on Lewis Research Center data only 



Frequency, f. Hz 


(c) = 0.8c^; V2=1.7 c^; 


= 1.0T ; 


T ^ = 3.8T . 
t,2 “ 
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